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1. Summary 

We present the first account of the structure -function relationships of a protein of 
the subfamily of urea/H + membrane transporters of fungi and plants, using 
Aspergillus nidulans UreA as a study model. Based on the crystal structures of 
the Vibrio parahaemolyticus sodium/ galactose symporter (vSGLT) and of the 
Nucleobase-Cation-Symport-1 benzylhydantoin transporter from Microbacterium 
liquefaciens (Mhpl), we constructed a three-dimensional model of UreA which, 
combined with site-directed and classical random mutagenesis, led to the identi- 
fication of amino acids important for UreA function. Our approach allowed us to 
suggest roles for these residues in the binding, recognition and translocation of 
urea, and in the sorting of UreA to the membrane. Residues W82, Y106, A110, 
T133, N275, D286, Y388, Y437 and S446, located in transmembrane helixes 2, 3, 
7 and 11, were found to be involved in the binding, recognition and/ or transloca- 
tion of urea and the sorting of UreA to the membrane. Y106, A110, T133 and Y437 
seem to play a role in substrate selectivity, while S446 is necessary for proper sort- 
ing of UreA to the membrane. Other amino acids identified by random classical 
mutagenesis (G99, R141, A163, G168 and P639) may be important for the basic 
transporter's structure, its proper folding or its correct traffic to the membrane. 



2. Introduction 

Urea is a readily available nitrogen source for bacteria, fungi and plants, occurring in 
nature as a product of animal, plant and microbial metabolism of nitrogenous com- 
pounds [1-6]. Additionally, together with ammonium nitrate it is one of the most 
widely used nitrogen-based fertilizers (http://espere.mpch-mainz.mpg.de/docu- 
ments/ACCENT/Edition08/texts%20material/english/fertilizer_usage_2001.pdf 
[7]). The internalization of this highly polar molecule necessitates specialized cellu- 
Electronic supplementary material is available lar transport systems. Fungi and plants possess a specific urea permease which is 
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and which constitutes a subfamily of H + -symporters within the 
sodium : solute symporter (SSS) superfamily. The latter is com- 
posed by hundreds of members, present in both eukaryotic 
and prokaryotic organisms, which usually function as Na + - 
cotransporters of a variety of substrates, such as sugars, vitamins, 
amino acids, inorganic ions, etc. (http://www.tcdb.org [8]). 

Members of this urea/H + symporters group have been 
functionally characterized in the fungi Saccharomyces cerevisiae 
(ScDUR3 [9]), Paxillus involutus (PiDUR3 [10]), Aspergillus 
nidulans (UreA [11]) and Candida albicans (CaDUR3 [12]), and 
in the plants Arabidopsis thaliana (AtDUR3 [13,14]) and Oryza 
sativa (OsDUR3 [15]). They share 35-74% of pairwise sequence 
identity and most of them include 14-15 predicted a-helical 
transmembrane segments (TMS). In A. nidulans, the urea trans- 
porter UreA has been identified as the only high-affinity urea 
transport system in this organism [11]. A second transport 
system may mediate facilitated diffusion at high urea concen- 
trations (more than or equal to 3mM) [11,16]. UreA is 
predicted to have 15 TMS [11]. Kinetic studies have established 
that it transports urea with a K m of 26 |xM, transport being 
inhibited by 2-thiourea (a urea toxic analogue) and acetamide. 
The transcription of the cognate gene, ureA, depends on the 
AreA GATA factor; it is repressible by ammonium and is not 
inducible by urea or its precursors. As for other permease 
genes, ureA expression is upregulated during the isotropic 
growth phase of the conidiospores, but, at variance with 
most permeases involved in the uptake of nitrogen sources, 
the gene is still under strict nitrogen catabolite repression in 
this phase of growth. A post-translational regulatory mechan- 
ism operates on UreA, which is endocytosed in response to 
ammonium. Putative orthologues of UreA are present in all 
sequenced Aspergilli. Moreover, three other paralogues of 
UreA have been identified in A. nidulans and are also present 
in other Aspergilli, though not all species genomes include 
all paralogues. These proteins may be low-affinity urea trans- 
porters, as ureA deleted strains (ureAA) only show residual 
growth when urea is used as sole nitrogen source. This fact 
has been attributed to a secondary energy-independent trans- 
port system, as mentioned above, or otherwise may be due to 
paralogues with different and unknown specificity, but 
which may residually transport urea [11]. 

The complex and hydrophobic nature of transporters 
makes their purification and crystallization difficult. However, 
three-dimensional modelling on solved crystal structures com- 
bined with random and directed mutagenesis can shed light on 
different aspects of a transporter's conformation and function, 
as illustrated clearly for the first time with homologues of the 
lactose permease LacY [17-22] and of the neuro transmitter- 
transporter prototype LeuT [23-26]. In the A. nidulans 
system, the power of this kind of approach has been illustrated 
by the thorough dissection of the high-affinity /high-capacity 
xanthine /uric acid transporter UapA [27,28]. 

The crystal structure of a member of the SSS superfamily, 
the sodium /galactose symporter of Vibrio parahaemolyticus, 
vSGLT, has been solved at approximately 3 A resolution 
[30,48]. It shows a five-helix inverted repeat motif (the 5-HIR- 
fold), also found in other bacterial transporters that do not 
share any sequence similarity and that transport a wide variety 
of substrates, like LeuT [31-33], Mhpl [29,34] and BetP [35,36] 
among others ([37,38] and references therein). A 'rocking- 
bundle' mechanism has been proposed for transporters 
bearing this structure, which implies a substrate-induced 
alternation between inward- and outward-facing hydrophilic 



cavities, and the opening and closing of specific gating amino 
acids [26,39,40]. We present below a structure -function analy- 
sis of UreA, based on comparative modelling, and random and 
site-directed mutagenesis. We identify mutations affecting sub- 
strate binding and/ or translocation and membrane sorting of 
the transporter. This constitutes the first such analysis of a 
transporter of this family. 

3. Results 

3.1. Rational design of site-directed mutagenesis 
and construction of mutant strains 

Following our previous characterization of UreA [11], we 
aimed to identify key residues involved in urea transport. We 
assumed that important residues for transport activity would 
be conserved in functionally characterized UreA orthologues 
(ScDUR3, CaDUR3, AtDUR3, PiDUR3 and OsDUR3) but not 
in its paralogue proteins (ANID_07373.1, ANID_02598.1 and 
ANID_07557.1), which do not significantly transport urea. 
Therefore, we aligned these eight homologous proteins against 
the primary sequence of UreA (figure 1). Residues fulfilling the 
above-mentioned criteria are Y106, A110, N275, A374, Y388, 
Y437 and S446. Following the predicted secondary structure 
of UreA [11], these residues are estimated to be located, res- 
pectively, in transmembrane helixes 3 (Y106, A110), 7 (N275), 
9 (A374) and 11 (Y437 and S446), while Y388 is located 
immediately after TMS9 (figure 1). 

We sought to construct a model of UreA comprising amino 
acids 64-525, spanning TMSs 2-13 (electronic supplemen- 
tary material, figure SI), based on the structure of vSGLT of 
V. parahaemolyticus [30]. This structure is representative of the 
transporter in a closed state with an inward-facing cavity. It 
must be noted that, owing to the low sequence identity 
between the crystallized protein and UreA (see Experimental 
procedures), the reliability of the model is limited to the posi- 
tioning of secondary structure elements and rough location 
of amino acids. The predicted TMSs of UreA modelled on 
vSGLT and Mhpl (see below) are shown in the electronic sup- 
plementary material, figure SI. TMSs 3, 7 and 11 are cognate to 
the segments where, in vSGLT, key residues of the substrate- 
binding site and translocation pathway are located. Thus, five 
out of seven conserved residues in characterized urea transpor- 
ters are located in these putatively important regions. In 
particular, Y106 and Y437 are predicted to face each other at 
the bottom of an inward-facing cavity (figure 2a, b), while 
A110 lines the cavity, located one helix-turn above Y106 
(figure 2a, c). N275 and S446 are located in helixes 7 and 11, 
which line the cavity, though neither of them faces it. The local- 
ization of residues Y106, A110, N275, S446 and Y437 supports 
their putative functional role, and thus they were selected for 
site-directed mutagenesis. Residues A374 and Y388, which 
are not predicted to localize to any of these helixes, were not 
included among the residues we chose to mutate. 

The residues discussed above were mutated as indicated 
in figure 1. In all cases, genomic DNA from a strain carrying 
a wild-type (wt) UreA-GFP fusion was used as the template, 
so that the subcellular localization of the resulting mutant - 
GFP fusions could be immediately assessed. Single copy 
transformants integrated to the ureA locus were obtained by 
transforming a ureAA strain, as described in Experimental 
procedures. All mutants were tested for utilization of urea 
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Figure 1. Multiple sequence alignment of UreA and homologues. Aligned sequences include A. nidulans UreA (Gl: 67516273), characterized orthologues in fungi and 
plants— ScDur3 of 5. cerevisiae (Gl: 51013791), AtDur3 of A thaliana (Gl: 9758728), PiDur3 of P. involutus (sequence kindly provided by Morel etal. [10]), CaDur3 of 
C. albicans (Gl: 68484979) and 0sDur3 of 0. sativa (Gl: 115483686)— and UreA paralogues AN7373 (Gl: 259483267), AN2598 (Gl: 259488035) and AN7557 (Gl: 
67901140). Putative TMSs and the intracellular helix between TMS3 and TMS4 (ICH3/4) are represented by grey and white rectangles above the sequence, respectively. 
For space reasons, only those TMSs, ICHs and loops including mutated residues are shown; omitted segments are represented by dashed lines. Fully conserved amino acids are 
shaded in black, and structurally conserved amino acids are shaded in grey. The mutations obtained in this work, whether by classical or directed mutagenesis, are boxed in 
black; residues affected by directed mutagenesis and conserved in UreA orthologues but not in its paralogues are denoted by a circumflex above the sequence, while those 
affecting amino acids conserved in all UreA homologues are denoted by an asterisk. The aligned residues affected by classical mutagenesis are denoted by a filled circle. 




Figure 2. Molecular model of UreA in the inward-facing, closed conformation, (a) Secondary structure elements are coloured by sequence number going continu- 
ously from red (amino terminal) to blue (carboxy terminal). The rough position of the membrane is indicated by the grey dotted lines with the extracellular surface 
on the upper part of the figure, (b) Same as (a), viewed from outside the membrane and rotated 90° along the membrane plane, (c) Solvent accessible surface 
calculated with a probe radius of 0.2 nm. The slice on the surface is taken perpendicularly to the membrane to show the inward-facing cavity. External 
and internal sides of the surface are yellow and grey, respectively. For clarity, only the semitransparent cartoon representation of helixes 3, 7 and 11 is 
shown, coloured as in (a). 



as nitrogen source, resistance to 2-thiourea and kinetics of 
urea uptake as described in Experimental procedures. Non- 
conservative changes yielded mutations Y106A, AllOD, 
N275I, Y437A and S446L. When grown on urea as sole nitro- 
gen source at concentrations of 0.625-5 mM, all substitutions 
described above result in phenotypes identical to those of the 
ureAA strain [11], showing residual growth similar to that 
seen in the absence of any nitrogen source, both at 37°C 
(figure 3a) and at 25°C (not shown). The resistance to 2- 
thiourea seen at concentrations of 0.625-5 mM is also com- 
parable to that observed for the ureAA strain except for 
Y106A (see below). Western blots show that this phenotype 
is not due to the absence of the protein (figure 3b). 

In accordance with growth tests, 14 C-urea initial apparent 
uptake rates for strains carrying mutations AllOD, N275I, 
S446L and Y437A were less than or equal to 2% of the 
wt (table 1). Mutation Y106A results in partial impairment 
of 14 C-urea uptake rates (approx. 14% of that of the wt). 
This is in agreement with growth tests on 2-thiourea, where 
sensitivity can be observed at concentrations more than or 
equal to 2.5 mM (figure 3a). Substrate binding was assessed 
by calculating K m values for urea and K t values for 2-thiourea, 
acetamide and guanidine. The latter is a structural analogue 
of urea that only binds UreA at concentrations more than 
or equal to 3 mM. Y106A increases the K m for urea and the 
Ki for 2-thiourea and acetamide by a factor of approximately 
2 (table 1), while the K{ for guanidine is diminished. It must 



be taken into account that phenotypic differences in growth 
plates can only be detected when K m/i values are higher 
than the substrate concentration used in growth tests, as is 
the case for Y106A. 

Conservative changes at the same selected positions were 
carried out to yield mutations Y106F, A110G, N275Q, Y437F 
and S446T (figure 1). When Y106 is changed to Phe, the 
mutant strain shows growth phenotypes on urea and 2-thiourea 
that are practically indistinguishable from the wt. 14 C-urea 
transport assays provide a K m for urea that is 1.6 times higher 
than those of the wt. The K { value for acetamide is slightly 
higher than that of the wt (312 + 17 jxM, against 237 + 14 |xM 
in the wt; table 1), while, as observed for Y106A, the K{ for gua- 
nidine is diminished. The initial urea uptake rate of mutant 
Y106F reaches only approximately 30% of that of the wt, 
which suggests that even if substrate affinity is only marginally 
affected, its transport across the membrane is further impaired. 
These differences are not detectable in growth tests. 

Mutations A110G, Y437F and N275Q affect similarly growth 
on urea and sensitivity to 2-thiourea. On urea, N275Q and Y437F 
are almost indistinguishable from the wt (figure 3d), whereas 
A110G shows a very subtle growth defect. The initial uptake 
rates of mutants Y437F and N275Q are approximately 50% 
and approximately 90%, respectively, as compared with the 
wt, while in the case of mutant A110G, the initial uptake rate 
attains approximately 40% of the wt. Urea K m found in strains 
carrying these mutations are 92 + 6 |jlM for A110G; 33 + 4 |xM 
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Figure 3. Characterization of strains bearing mutations in residues conserved in functionally characterized UreA orthologues. (a) Growth phenotypes of mutant UreA 
strains at 37°C on urea as nitrogen source or on 2-thiourea with NaN0 3 10 mM as nitrogen source to test resistance to the analogue. Growth on ammonium 5 mM 
and nitrate 10 mM are used as controls. A wt and a ureAA strain are shown as positive and negative controls, respectively. Similar results were obtained at 25°C 
(not shown), (b) Western blot analysis of total protein extracts of UreA-GFP mutants probed with anti-GFP antibody. Cultures were grown for 14-16 h at 25°C in 
derepressed conditions (proline as sole nitrogen source). The low mobility band corresponds to intact UreA-GFP and the high mobility band to free GFP (see text). 
Antibody against actin was used as an internal control of loading, (c) Epifluorescence and confocal microscopy of mutant bearing the S446L mutation grown in 
derepressing conditions, showing the retention of the UreA-GFP fusion in the ER, seen as perinuclear rings (indicated by arrows). Upper panels show the local- 
ization of the GFP signal. Wild-type UreA-GFP localization is shown as control. Middle panels show the localization to the nucleus of a histone H1 (HhoA)-mRFP 
fusion. Lower panels display the colocalization of GFP and RFP signals. Scale bar, 10 |xm. 



Table 1. Summary of UreA mutations obtained by site-directed mutagenesis and their effect on subcellular localization and functionality. TMS, transmembrane 
segment; ICH3/4, intracellular helix between TMS3 and TMS4; M, localization to the membrane; ER, retention in the endoplasmic reticulum; n.a., not applicable. 

+++, growth phenotype on urea as that of the wt; , growth phenotype on urea as that of the ureAk; ++-, + — , intermediate growth 

phenotypes between that of the wt and that of the ureAk strain. V (%): initial uptake rates, expressed as % of the wt in pmoles min~ 1 /10 7 conidiospores. 
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According to the constructed three-dimensional models. 



for N275Q and 88 + 9 |jlM for Y437F (table 1), compared with a 
value of 26 + 2 |xM for the wt, which agrees with the observed 
changes in initial uptake rates. On 2-thiourea, A110G and 
Y437F show a moderate sensitivity at concentrations of 5 mM, 
while N275Q is already sensitive at concentrations of 2.5 mM 
(figure 3a). All three mutations result in K { values for 2-thiourea 
that are approximately twice as high as those of the wt and 
exceed the lowest concentrations used in plate tests (625 |xM; 
not shown). As mentioned above, this may explain why pheno- 
typic differences became visible for these mutants. Additionally, 
Ki values suggest that Y437F results in an augmented apparent 
affinity for acetamide (102 + 15 |xM compared with 237 + 
14 fxM in the wt; table 1), while A110G and N275Q result in 
reduced affinity for acetamide (495 + 29 |jlM and 352 + 
26 |xM, respectively, compared with 237 + 14 |xM in the wt; 
table 1). S446T does not affect the phenotype in growth tests 
nor does it change significantly the uptake kinetic constants. 

Epifluorescence microscopy was used to investigate the sub- 
cellular localization of GFP-tagged UreA versions bearing all the 
above-mentioned mutations. With the exception of S446L, none 
of these mutations (Y106A/F, A110D/G, N275I/Q, Y437A/For 
S446T) affect the membrane localization of UreA (electronic 
supplementary material, figure S3). In non-repressive condi- 
tions (proline as nitrogen source), the GFP signal is seen 
predominantly in the plasma membrane and the septae, and 



to a lesser extent in cytoplasmic organelles previously identified 
as endosomes and vacuoles [11]. S446L results in a completely 
different localization, with fluorescence accumulating in ring- 
shaped structures around the nucleus (figure 3c), a pattern 
typical of endoplasmic reticulum (ER)-resident proteins 
[41,42]. Fluorescence also accumulates in vacuoles. The peri- 
nuclear localization was confirmed with a strain carrying 
UreAS446L-GFP and a Hl-mRFP fusion [43,44] (figure 3c). 
This localization suggests that changing S446 into an L may be 
detrimental for proper folding of UreA or for its trafficking 
towards the plasma membrane, causing retention in the ER. 
Western blots show a diminished amount of total UreA protein 
as compared with the wt (intact UreA-GFP + free GFP, the 
latter resulting from UreA-GFP turnover; figure 3b). 

3.2. Directed mutagenesis of residues in the putative 
substrate pathway 

We then proceeded to analyse in more detail the conserva- 
tion of residues lining the putative pathway of urea through 
the transporter. We reasoned that residues involved in this 
function should (i) be conserved among the family of urea 
transporters, (ii) be solvent exposed in the structural model, 
and (iii) possess polar physicochemical properties to interact 
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Figure 4. Analysis of strains bearing mutations in the putative intramolecular urea pathway, (a) Location of the mutated residues, (i) A perspective view into the 
inward-facing cavity. W82, T133, Y388 and D286 are shown from the cytoplasmic side. Y106 and Y437, which close the cavity, are shown for reference. The colouring 
of the cartoon is the same as in figure 2. The additional intracellular helix where T133 is located is shown as semitransparent in orange, (ii) A side view of the 
transporter showing the W82, T133, Y388 and D286 residues in the predicted urea path. For clarity, we show only the helixes where the residues chosen for directed 
mutagenesis are located, (b) Growth tests of strains bearing mutations in W82, T133, Y388 and D286 residues. Growth conditions were as described in figure la. 
yA + and yA2 wt strains were used as positive controls so that each mutant strain could be compared with the relevant one. A ureAA strain is shown as a negative 
control, (c) Western blot analysis on total protein extracts of UreA-GFP mutants probed with anti-GFP antibody. See legend of figure lb for details. 
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with the urea molecule. This analysis led to the identification 
of four amino acids: W82, T133, D286 and Y388 (figure 4a). 
W82 is situated in TMS2, T133 in the intracellular helix 
found between TMS3 and TMS4 (ICH3/4), D286 in TMS7, 
and Y388 at the end of TMS9 (electronic supplementary 
material, figure SI). Interestingly, the latter is among the 
amino acids identified as putatively involved in transporter 
functionality because of its conservation in functionally charac- 
terized UreA orthologues but not in its paralogues (see above). 

These four residues were mutated to alanine and con- 
servative changes were also introduced at these positions. 
Phenotypic analysis of mutant strains, apparent urea transport 



rates and uptake competition assays were carried out as men- 
tioned above. In agreement with structural predictions, W82A, 
T133A, D286A and Y388A result in loss of function of UreA, 
because mutants carrying them are unable to grow on urea at 
concentrations of 0.625-2.5 mM and are resistant to 2-thiourea 
when present at the same concentration range (figure 4b). 
Accordingly, initial urea uptake rates are in all cases almost 
undetectable (less than 2%; table 1). 

Mutants carrying the D286E and Y388F conservative 
changes grow normally on urea as sole nitrogen source and 
show normal sensitivity to 2-thiourea (figure 4b). Strains carrying 
these mutations show 100% apparent urea uptake rates and K m /i 




G99S R141C A163P 



Figure 5. Analysis of mutants isolated by classical random mutagenesis, (a) Structural context of mutants isolated by classical random mutagenesis, (i) The UreA 
molecule is shown as helix-succession representation, coloured as in figure 2. The residues affected by the mutations are shown with a space-filling representation. 
For reference, Y106, A110, N275, Y437 and S446 are shown by stick representation, (ii) Same as (i), but seen from the extracellular side. Stick representations of Y106 
and Y437 are included for reference, (b) Epifluorescence microscopy in greyscale inverted mode, showing in vivo subcellular expression of mutant UreA-GFP fusions. 
Arrows signal perinuclear ER membrane rings, (c) Western blot analysis on total protein extracts of UreA-GFP mutants probed with anti-GFP antibody. See legend of 
figure 3b for details. 



affinity constants for urea, 2-thiourea, acetamide and guanidine 
that are not significantly different from those of the wt strain 
(table 1). T133S is phenotypically indistinguishable from the wt 
(figure 4b). Nevertheless, this mutation results in an initial urea 
uptake rate of approximately 70% of that of the wt. This may 
be explained by a reduced affinity for urea (K m 80% higher com- 
pared with the wt). K { for 2-thiourea and acetamide is also 
increased (approximately 30%) for these two substrates (table 
1). Finally, mutant W82F behaves phenotypically on urea as a 
partial loss-of-function mutant. By contrast, no phenotypic 
differences are obvious on 2-thiourea compared with the wt 
(figure 4b). The apparent urea uptake capacity for this mutant 
is approximately 50% of that of the wt. However, its K m for 
urea is not significantly different to that of the wt (table 1), 
which means that even if the substrate could bind efficiently, 
the mechanism of translocation of urea across the membrane is 
somehow hindered by the W82F mutation. Moreover, K { 
values for other substrates are slightly increased, which points 
to a moderate role of this residue in substrate selectivity. Epifluor- 
escence observations of this set of mutants show that in all of 
them UreA localizes normally to the membrane (electronic sup- 
plementary material, figure S3). The cognate mutations do not 
produce drastic differences in protein concentrations in total 
preparations (figure 4c), as demonstrated by Western blots. 

3.3. Mutations obtained by random mutagenesis 

As a complementary strategy, the isolation of ureA mutants by 
classical random mutagenesis (as described in Experimental 
procedures) and the characterization of these and previously 
isolated spontaneous mutants were carried out. Mutated 
residues are summarized in table 2 and shown in figure 5a. 



Table 2. Summary of UreA mutations obtained by classical random 
mutagenesis. TMS, transmembrane segment; TMS2/TMS3, loop between 
TMS2 and TMS3; ICH3/4, intracellular helix between TMS3 and TMS4; C-ter, 
C-terminal region; M, localization to the membrane; ER, retention in the 
endoplasmic reticulum; n.d., not determined. 
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Mutation G168D [11,16] is predicted by TMHMM [45] 
and our structural model to affect the last residue in TMS4, 
while mutation P639R [11] affects the first residue in the 
C-terminal region of the protein (figures 1 and 5a). This 
region is not included in our structural model. Mutants 
bearing G168D or P639R show only a residual growth on 
urea, similar to that of the ureA A strain [11]. 14 C-urea trans- 
port assays agree with these results [11]. Epifluorescence 
microscopy observations of GFP fusions of these mutant pro- 
teins show the characteristic localization of ER-retained 



proteins, in ring-shaped structures around nuclei (figure 5b), 
whereas in the case of mutation G168D some fluorescence is 
also seen in the plasma membrane. 

Five novel mutations were selected for further work, as 
described in Experimental procedures: G99S (predicted to be 
in the extracellular loop between TMS2 and TMS3), R141C 
(ICH3/4), A163P (in TMS4), W480 to a STOP codon (end of 
TMS12) and Q570 to a STOP codon (predicted to be in the extra- 
cellular loop between TMS13 and TMS14) (figures 1 and 5a). 
Mutants bearing any of these changes behave like the ureAA 
mutant in terms of growth phenotypes on urea and 2-thiourea, 
when incubated at 37°C (electronic supplementary material, 
figure S4) or 25 °C (not shown). No differences were observed 
in terms of initial uptake rates with respect to the ureAA strain 
(not shown). 

Epifluorescence microscopy showed that in the A163P 
mutant, UreA localization is indistinguishable from that of a 
wt strain, staining mainly the membrane and septae, while 
mutations G99S and R141C result in UreA retention in the 
ER (figure 5b). Fluorescence is also present in vacuoles in 
both mutants, a fact that could imply an ER-to-vacuole degra- 
dation route [46] or transient plasma membrane localization 
followed by degradation due to protein instability. 

Western blots of total protein extracts (figure 5c) show that 
neither diminution of intact UreA-GFP levels nor any increase 
in degradation (as shown by the relative level of free GFP) 
occurs in these mutants. The increase in protein level observed 
for G168D could be due to an increase in the intrinsic stability 
of the mutant protein or to reduced vacuolar turnover due to 
ER retention [47]. 

3.4. Modelling of UreA in the outward-facing cavity 
conformation 

Results obtained for strains bearing mutations Y106F and W82F 
show that there is a reduction in the initial transport rate which 
does not correlate with K m and K t changes (table 1). This may 
suggest that urea can bind to the protein but the mutation some- 
how impairs the conformational rearrangements of the protein 
along the translocation process. To acquire a different structural 
perspective on the role of these residues, we undertook the con- 
struction of a model of UreA based on the structure of the 
Nucleobase-Cation-Symport-1 benzylhydantoin transporter 
from M. liquefaciens (Mhpl; PDB code 2JLN [29]). Although 
Mhpl shows less amino acid sequence identity with UreA 
than vSGLT (see Experimental procedures), this model was con- 
structed to obtain an insight into the possible role of the amino 
acids identified in the outward closed conformation. As it can 
be expected, the models constructed using structural templates 
in inward closed and outward closed conformations are signifi- 
cantly different. Indeed, the positions and interactions of most of 
the residues undergo substantial changes upon the conformatio- 
nal transition from the inward- to the outward-facing cavity 
(figure 6; electronic supplementary material, movie). 

The large rearrangements undergone by Y106 and W82 
comparing inward- and outward-facing conformations lead 
us to speculate that changing any of them into Phe may 
impair the conformational transition, without substantially 
affecting the binding of urea in the outward facing cavity. The 
different perspective provided by this structural model may 
also suggest a possible role for R141. In the conformation pre- 
senting the outward-facing cavity, R141 is in good position to 



form a salt bridge interaction with E123 (electronic supplemen- 
tary material, figure S5), which is conserved only in AtDUR3 
and OsDUR3 (figure 1). This salt bridge would be disrupted 
by the R141C mutation, causing a structural modification that 
may ultimately lead to the retention of UreA in the ER. 

4. Discussion 

In this work we present, based on a strategy combining 
site-directed and classical random mutagenesis with three- 
dimensional modelling, the first structure -function study of a 
protein of the group of urea/H + transporters of fungi and 
plants, the A. nidulans UreA. Although the reliability of the mod- 
elling is limited by the relatively low amino acid identity with 
structurally characterized homologues, the models obtained 
allowed pinpointing specific amino acids, which, together with 
classical random mutagenesis, led to the identification of specific 
residues important for the function of UreA, some of which are 
crucially involved urea uptake by UreA. 

We selected residues W82, Y106, A110, T133, N275, D286, 
Y388, Y437 and S446 for site-directed mutagenesis. Except for 
N275 and S446, these residues are predicted by both our 
three-dimensional models to be exposed to the cavity and 
thus would be able to establish contacts with the substrate. 
N275 and S446 are predicted to be in helixes that line the 
cavity, but are not exposed to it. 

Interestingly, our model suggests that W82, Y106, A110, 
N275, D286, Y437 and S446 are located in TMSs 2, 3, 7 and 
11. These helixes correspond to helixes 2, 3, 7 and 11 in 
vSGLT, where they are involved in substrate binding and 
recognition [30,48]. 

Our data indicate that the above-mentioned residues would 
be involved in urea transport at different steps of the substrate 
translocation process. With the exception of S446, non- 
conservative mutations of these residues do not alter the 
sorting of UreA mutant proteins towards the membrane; 
transport impairment must therefore be due to reduced 
binding or translocation of the substrate across the membrane. 
With the exception of D286E, Y388F and S446T, all conservative 
changes introduced in UreA affect urea transport. A reduction 
in initial transport rates correlates well with changes in K m 
and K ir except for strains carrying mutations W82F and 
Y106F, whose K m/i values do not fully explain the impair- 
ment in uptake rates. It can then be speculated that these 
mutations must also affect other aspects of the transporter's 
function. Finally, changes in K m and K { for substrates and 
inhibitors resulting from mutations Y106F, A110G, T133S 
and Y437F suggest involvement of these residues in UreA 
substrate selectivity. 

5-HIR transporters are characterized by the presence of 
filters or gates that determine the open or occluded states in 
the inward- and outward-facing conformations, which may 
have a role in substrate selectivity [37,40]. The location of 
residues Y106 and Y437, which, according to the model, are 
predicted to face each other in the inward closed conforma- 
tion, could suggest a role of these two aromatic residues as a 
filter gate [37,49,50]. Moreover, these two residues lay in a 
position that almost superimposes with residues shown to be 
part of the outside gate in vSGLT. The fact that these two 
residues seem to play a role in substrate selectivity reinforces 
this idea. Our molecular model of the UreA in the inward- 
facing, closed conformation suggests that selective permeation 
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Figure 6. Molecular model of UreA in the outward-facing conformation, (a) Secondary structure elements and amino acids are coloured as in figure 2. The approxi- 
mate position of the membrane is indicated by the grey dotted lines with the extracellular side on the upper part of the figure, (b) Same as in (a), rotated 90° 
along the membrane plane and seen from the extracellular side, (c) Solvent accessible surface calculated with a probe radius of 0.2 nm. The slice on the surface is 
taken perpendicular to the membrane to show the outward-facing cavity. External and internal sides of the surface are yellow and grey, respectively. For clarity, we 
show only the helixes where the residues chosen for directed mutagenesis are located, (d) Perspective view into the outward-facing cavity from the intracellular side 
(compare with figure 5a). The colour of the cartoon is the same as figure 2, and the intracellular helix where T133 is located is shown as semitransparent. 



may be regulated in a similar fashion to that deduced from 
the crystal structure of a bacterial channel-like urea trans- 
porter, dvUT [51]. In this case, urea would pass through a 
selectivity filter where a couple of facing Phe residues 'sand- 
wich 7 the substrate, so that the aromatic rings stabilize the 
partial positive charge on the amide nitrogen by amide -PI 
stacking interactions. 

Mutation S446L is unique among the site-directed mutations 
studied in this work, in the sense that it is the only one among the 
rationally designed mutations that abolishes the correct traffick- 
ing or insertion of the mutant protein into the plasma membrane, 
while the substitution by threonine maintains normal function. 
This suggests that a non-conservative change may lead to an 
aberrant folding of the mutant protein. 

Classical random mutagenesis led to the isolation of 
mutants in TMS4, in the loops between TMS2 and TMS3, 
TMS3 and TMS4 and TMS4, and TMS5, and also at the begin- 
ning of the predicted C-terminal region of the protein. Except 
for mutation A163P, all other changes lead to the retention of 
the mutant protein in the ER. This suggests a relevant role of 
these residues in the establishment of proper folding of the 
protein or in its correct traffic to the membrane. It is worth 
noting that in support of a role in the determination of the 
basic transporter's structure, all the residues identified by 
classical random mutagenesis show a high degree of conser- 
vation not only among UreA orthologues but also, for a 
number of them, among SSS members. Residues G99 and 



R141 are conserved in UreA orthologues and in all members 
of the SSS family characterized so far (figure 1; electronic 
supplementary material, figure S2), while G168 is partially 
conserved in characterized SSS members, where a serine is 
found in the same position in S. cerevisiae DUR3 and a proline 
in NIS (human sodium-iodide symporter; electronic sup- 
plementary material, figure S2). P639 is conserved in all 
UreA orthologues, but not in SSS members. One possible 
hypothesis is that the mutation of this residue abolishes a 
kink at the carboxy terminus of TMS15, which may be impor- 
tant for the correct orientation of the C-terminal region of 
UreA. C- and N-terminal tails are characteristic of eukaryotic 
transporters bearing 5-HIR structures [27,28,47,52], while in 
bacterial transporters of this family these regions are not pre- 
sent. A163 is the less conserved residue identified by random 
mutagenesis (figures 1; electronic supplementary material, 
S2); it is substituted by a serine in ScDUR3, PiDur3 and 
CaDur3, and is not conserved at all in non-UreA homologues 
of the SSS. It could be speculated that a proline substitution 
would introduce a kink in helix TMS4. 

The work presented here underlines the opportunity to 
combine modelling with site-directed and random mutagen- 
esis analysis. While directed mutagenesis suggested the 
important role of residues for substrate binding and translo- 
cation, random mutagenesis identified residues involved in 
protein translocation to the membrane and the possible role 
of a salt bridge to stabilize the conformation of the protein. 
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Finally, our work could not provide, in the absence of other 
specific information, any insight on the amino acids or regions 
of UreA involved in proton symport. In previous work [11], we 
performed inhibition experiments in the presence of the 
H + -uncoupler CCCP and the H + -ATPase inhibitor DCCD, 
which completely abolish urea uptake by UreA. These exper- 
iments, despite being indirect, show that, like most fungal 
transporters, UreA operates as an H + /symporter. However, 
there is no way of distinguishing mutations resulting in an 
impairment of H + / symport from other loss-of-f unction 
mutations, so no conclusions can be drawn from inhibition 
experiments of mutants. In this work, all mutants that con- 
served some activity, when tested in the presence of CCCP, 
showed abolished uptake similarly to the wt, therefore these 
results were not shown. 

5. Experimental procedures 

5.1. Strains, media and transformation procedures 

Standard complete and minimal media (MM) for A. nidulans 
were employed [53,54] (http://www.fgsc.net). Supplements 
were added when necessary at standard concentrations 
(http:// www.gla.ac.uk/ ibis/ molgen/ aspergillus/ supplement, 
html). Aspergillus nidulans strains used in this study are listed 
in the electronic supplementary material, table SI. Gene sym- 
bols are defined in http://www.gla.ac.uk/ibls/molgen/ 
aspergillus/loci.html. Urea (0.625-5 mM), NaN0 3 (10 mM), 
ammonium L(+)-tartrate (5-10 mM) or proline (5 mM) were 
used as sole N-sources. 2-thiourea was used in concentrations 
of 0.625-5 mM. Aspergillus nidulans transformation was carried 
out as in [55]. 

5.2. Generation of ureA mutants 

5.2.1. Site-directed mutagenesis 

Site-directed mutagenesis on the ure A:\gfp gene was per- 
formed by the Fusion-PCR technique [55] using KAPA HiFi 
DNA polymerase (KAPA Biosystems) with primers Ure5-F, 
Ure3-R and complementary ones carrying the desired substi- 
tution (electronic supplementary material, table S2). DNA 
extracted from a ureA::gfp::AFpyrG strain (MVD 10A) was 
used as template. A resulting 7 kb fusion product was ampli- 
fied with nested primers Ure5-N and Ure3-N and purified 
with the Gene JET Gel Extraction Kit (Thermo Scientific). 
The resulting construct was transformed in a ureAk::riboB 
pyrG89 pyroA4 riboBl nkuAA::argB veAl strain (MVD 13A). 
Protoplasts were plated on selective (containing riboflavin, 
but not uridine and uracil) 1 M sucrose MM and incubated at 
37°C. Transformants were purified on selective MM and its 
ability to grow both in different concentrations of urea and 
2-thiourea was tested. Southern blots were performed in 
order to check the integration to the ureA locus and to assess 
the number of construction copies integrated in the genome. 
Probes were labelled with the DIG DNA labelling and detec- 
tion kit (Roche Applied Science). Sequencing of the resulting 
mutants was performed at Macrogen, Inc. (Seoul, Korea). 

5.2.2. Classical random mutagenesis 

Classical random mutagenesis was performed by incubating a 
10 ml spore suspension from strain MVD 10 A with 20 julI of 



4-nitroquinoline 1-oxide (1 |xg ml -1 ) for 1 h at 37°C. Selection 
of 2-thiourea-resistant strains was performed on medium 
containing 1.25 mM 2-thiourea using NaN0 3 (10 mM) as a 
nitrogen source. Five mutants were selected for further 
analysis. Crossover with an ureAl strain (MVD 100) was per- 
formed in order to verify that the produced mutation was in 
the ureA locus. Sequencing of ureA mutants was performed 
at Macrogen, Inc. 

5.3. Radiolabeled urea uptake measurements 

[ 14 C]-urea uptake in MM was assayed in germinating conidio- 
spores of A. nidulans at 37°C, pH 6.8, as previously described 
[11]. Initial velocities were measured at a time point of incu- 
bation where velocity increases linearly (2 min), using a 
concentration of 1.0 |xM [ 14 C]-urea, and are expressed as % of 
the wt in pmoles min _1 /10 7 conidiospores. K m/i values were 
obtained directly by performing and analysing uptakes 
(Prism 3.02; Graph Pad Software, Inc.), using labelled urea at 
0.5-50 |xM, or in the presence of 1.0 |xM [ 14 C]-urea and of var- 
ious concentrations (0.5-3000 |jlM) of non-labelled substrates. 
K { values were calculated from the Cheng and Prusoff equation 
Ki ■ = IC 50 /(1 + (L/K m ), in which L is the permeant concen- 
tration. Transport assays were carried out in at least three 
independent experiments, with three replicates for each con- 
centration or time point. Background uptake values were 
corrected by subtracting values measured in the deleted 
mutants (UreAA). Standard deviation was less than 20%. 
Radiolabeled [ 14 C]-urea (55.0 mCi mmol" 1 ) was purchased 
from Moravek Biochemicals, Brea, CA, USA. 

5.4. Protein extraction and Western blot 

Total protein of 200 mg of grinded mycelia was extracted as 
described by Apostolaki et ah [56] and protein concentration 
was determined by the Bradford assay. Total proteins (50 |xg) 
were separated by SDS-PAGE (10% (w/v) polyacrylamide 
gel) and electrob lotted (OmniPAGE Electrob lotting Units, 
Cleaver Scientific Ltd) onto a nitrocellulose membrane 
(0.45 |xm, Thermo Scientific). Membrane was then incubated 
in stripping buffer (10 mM Tris-HCl pH 6.8, 2% SDS and 
100 mM p-mercaptoethanol) for 15 min at 55°C (as described 
by Kaur & Bachhawat [57]). Membrane was then treated with 
5% non-fat dry milk, and immunodetection was done with 
a primary mouse anti-GFP monoclonal antibody (Roche 
Applied Science), or a mouse anti-actin monoclonal (C4) anti- 
body (MP Biomedicals) and a secondary sheep anti-mouse 
IgG HRP-linked antibody (GE Healthcare). Blots were devel- 
oped using the Amersham ECL Western blotting detection 
reagents and analysis system (GE Healthcare), and images 
were acquired using the GeneSys software from the GBoxChemi 
XT4 System (Syngene). 

5.5. Epifluorescence microscopy 

Samples for fluorescence microscopy were prepared as 
described previously [58]. In brief, samples were incubated 
directly on sterile coverslips protected from light in liquid MM 
with proline (5 mM) as nitrogen source and appropriate 
supplements, at 25°C for 14-16 h. For epifluorescence micro- 
scopy, samples were observed in an Olympus inverted 
microscope CKX31 with a U-MNIBA3 (for GFP) or U-MWIG3 
(for RFP) filter, photographed with a Hamamatsu Orca Er 



camera and processed using Image Pro v. 6.0 software. For 
confocal microscopy, samples were observed in a Leica inverted 
microscope DMI6000CS with an 13 (for GFP) or N2.1 (for RFP) 
filter. Images were acquired and processed using LASAF 
v. 2.20 software. The resulting images were processed by 
Adobe Photoshop software. Microscopy facilities belonging to 
the Cellular Biology Platform, Institut Pasteur de Montevideo, 
were used. 

5.6. Bioinformatics and modelling of UreA 

Sequences were obtained from the Aspergillus Comparative 
Database, Broad Institute (http://www.broad.mit.edu/ 
annotation/ fungi /aspergillus), the Saccharomyces genome 
database (http://www.yeastgenome.org), the Candida 
Genome Database (http://www.candidagenome.org), The 
Arabidopsis Information Resource (http://www.arabidop 
sis.org), the Oryza sativa genome database (http://www. 
plantgdb.org/ OsGDB) and from the National Center for Bio- 
technology Information (http: / / www.ncbi.nlm.nih.gov) . 
PiDUR3 sequence was kindly provided by Morel et al. [10]. 
Multiple sequence alignments were carried out in MEGA v. 
5.0 software [59] with Muscle [60]. Homologues of UreA 
with known structure were identified by performing a 
remote homology search on the Protein Data Bank (http:// 
www.pdb.org) using Hidden Markov model profiles with 
Hhpred [61]. The best structural template in the closed 
inward-facing conformation corresponds to the structure of 



vSGLT of V. parahaemolyticus [30] with a sequence identity 
of 18%, and E- and p-values of 3.77 x 10" 65 and 1.2 x 
10 ~ 69 , respectively. The best structural template in the 
closed outward-facing conformation corresponds to the 
structure of Mhpl from M. liquefaciens [29] with a sequence 
identity of 14%, and E- and pvalues of 5.6 x 10~ 09 and 
1.8 x 10 ~ 13 , respectively. Molecular models on both confor- 
mations comprising amino acids 64-525, spanning TMS 
2-13, were built using those templates and the sequence 
alignments obtained from HHpred with MODELLER 9.12 
[62]. In both cases, 1000 models were generated, from 
which the best 10 were selected on the base of the Modeller's 
Objective Function. The stereochemical properties of the 
remaining models were assessed using Procheck [63]. The 
models with highest number of amino acids in most favou- 
red regions and higher overall G-factors were selected as 
representatives of the inward- and outward-facing confor- 
mations (electronic supplementary material, table S3). In 
both cases, they are comparable or better than expected for 
medium-quality X-ray structure. 
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